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Resumen 


Los cálculos oportunos de la evapotranspiración de referencia diaria (ETo) 
cuantifican y planifican el uso de la disponibilidad hídrica y, en 
consecuencia, determinan un consumo diario, más en las condiciones 
donde la falta de información climática interrumpe el rendimiento de los 
diversos modelos para estimar la evapotranspiración de referencia diaria 
en una determinada región. El objetivo de este trabajo es evaluar el 
comportamiento de las ecuaciones de Hargreaves-Samani, Priestley y 
Taylor, García y López, y Linacre en el cálculo de la evapotranspiración de 
referencia diaria (ETo) respecto a la ecuación de FAO Penman-Monteith 
(PM), utilizando la información diaria de cinco estaciones próximas a la 
zona productora de arroz del Valle del Cauca, Colombia. Se estimó la 
evapotranspiración de referencia diaria de 28 años consecutivos; se 
utilizaron los indicadores estadísticos: coeficiente de determinación (R2); 
la raíz cuadrada del cuadrado medio del error (RMSE); eficiencia de Nash- 
Sutcliffe (NSE), y diagramas de cajas para definir la bondad de los 
modelos utilizados. Los modelos Priestley-Taylor (PT) y Linacre (LN) han 
exhibido los resultados más afines al modelo estándar. Las estaciones 
Ing. Bengala y Ciat Quilichao se exhiben como las áreas que mejores 
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comportamientos poseen de todos los modelos con el índice del 
coeficiente de determinación (R?); todas las estaciones presentaron 
resultados aceptables con la raíz cuadrada del cuadrado medio del error 
(RMSE) para las diversas comparaciones realizadas, no así con la 
eficiencia de Nash-Sutcliffe. Existe evidencia limitada del uso diario del 
modelo García y López, siendo este método y el de Hargreaves-Samani 
los menos ideales para la estación Independencia. 


Palabras clave: modelos de evapotranspiración, arroz, datos diarios, 


evapotranspiración de referencia. 


Abstract 


The opportune calculations of the daily reference evapotranspiration 
(ETo) quantify and schedule the use of the water availability and 
consequently determine a daily consumption, mostly in the conditions 
where the lack of climatic information interrupts the performance of 
diverse models to estimate the daily reference evapotranspiration in a 
certain region. Therefore, the objective of this work is to evaluate the 
behavior of the equations of Hargreaves-Samani, Priestley and Taylor, 
García and López, and Linacre in the calculation of daily reference 
evapotranspiration (ETo) for the equation of FAO Penman-Monteith (PM) 
using daily information from five stations nearthe rice production area of 
Valle del Cauca, Colombia. The daily evapotranspiration of 28 consecutive 
years was estimated, the results were confirmed by the statistical 


indicators: Determination Coefficient (R?); the square root of the mean 
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square of the error (RMSE); Nash-Sutcliffe Efficiency (NSE), and box 
diagrams. The Priestley-Taylor (PT) and Linacre (LN) models have shown 
the closest results to the standard model. Bengala and Ciat Quilichao 
stations are shown as the areas with the best behavior of all the models 
with the coefficient of determination index (R?); all the stations presented 
acceptable results with the square root of the mean square of the error 
(RMSE) for the different comparisons made but not with the Nash-Sutcliffe 
Efficiency (Nash-Sutcliffe Efficiency). However, there is limited evidence 
of the daily use of the Garcia and Lopez model, being this method and 
that of the Hargreaves-Samani the least ideal for the Independence 


station. 


Keywords: Evapotranspiration models, rice, daily data, reference 


evapotranspiration. 
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Introducción 
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Las estimaciones de evapotranspiración de referencia (ETo) han sido 
ampliamente estudiadas y utilizadas en los fundamentos de ingeniería de 
riego para definirlas necesidades de agua en los cultivos. De acuerdo con 
Sánchez (2000), el término de evapotranspiración de referencia (ETo) ha 
sido usado para describir los procesos totales de trasferencia de agua a 
la atmósfera desde una vegetación con características bien definidas. En 
la agricultura, la ETo es una de las variables elementales para promover 
un uso más eficiente del agua, por esto la necesidad de tener 
estimaciones adecuadas (Back, 2008), ya que juegan un papel esencial 
en el balance hídrico para el manejo de cuencas y para la estimación de 
los requerimientos de riego, la programación del riego y estudios sobre el 


cambio climático. 


El sur del Valle del Cauca-Colombia posee dos épocas de siembra a 
nivel secano con 3 750 ha productoras de arroz (Dane € Fedearroz, 
2017); la cuantificación de la ETo diaria para la zona a nivel diario es 
necesaria para mejorar el manejo del cultivo, conocer sus necesidades 
hídricas y realizar un uso eficiente del agua (Martínez € Pérez, 2006); 
pero para estimar dicha variable se carecen de registros diarios para 
utilizar el modelo Penman-Monteith (PM) propuesto por Allen, Pereira, 
Raes y Smith (1998). En Colombia existen algunos trabajos relacionados 
con la estimación de la ETo a diferentes escalas: Barco et al. (2000) 
concluyeron que Penman y Morton son métodos difíciles de implementar 
debido a la información climática requerida; la ecuación Cenicafe 
desarrollada por Jaramillo (2006) permite calcular la evapotranspiración 
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(ETo) de forma sencilla, dicha ecuación solo depende de la elevación sobre 
el nivel del mar; la ETo estimada para todo el país por Poveda et al. (2007) 
indica que las ecuaciones presentadas por Turc y Morton son las más 


convenientes para estas evaluaciones. 


La tendencia de la estimación de este tipo de variable meteorológica 
ha indicado que se reportan más resultados obtenidos a escala anual y 
mensual; no obstante, se han identificado en Scopus algunos reportes 
significativos sobre la estimación de la evapotranspiración de referencia 
diaria utilizada en la estimación de la ET del cultivo del arroz publicados 
en 2017 en Taiwán y en 2019 en Senegal. Al recurrir a los ejemplos, Yao 
et al. (2017) obtuvieron los mapas de la ETo promedio de 10 días a partir 
de imágenes satelitales meteorológicas, comparada con la ETo derivada 
del modelo PM de la FAO (Food and Agriculture Organization of the United 
Nations) (Allen et a/., 1998); Djaman, Sall, Sow, Manneh e Irmak (2019) 
determinaron los datos climáticos medidos en el área del arroz de riego 
que se utilizan en la estimación de la ETo mediante el modelo Penman- 
Monteith en condiciones en las que no existen registros de las estaciones 


meteorológicas. 


La elección de un modelo o método de cálculo depende de la 
disponibilidad de los datos meteorológicos, así como de la exactitud y 
precisión del modelo para una región dada (De-Carvalho et a/., 2013). 
Por ende, La FAO ha adoptado el modelo de Penman-Monteith 
principalmente por dos razones: a) se trata de un modelo físico aplicable 
no sólo a una zona en específico, y b) los datos que se usan se pueden 
obtener de estaciones meteorológicas, observatorios o mediciones a 
través de lisímetros (López, Hess, € White, 2009). Expuesto todo lo 
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anterior, existe para el valle del Cauca por lo menos un modelo que estime 
con un error bajo la ETo a nivel diario con respecto a un modelo estándar. 
El objetivo de este trabajo fue evaluar el comportamiento de las 
ecuaciones de Hargreaves-Samani, Priestley-Taylor, García y López, y 
Linacre en el cálculo de la evapotranspiración de referencia diaria (ETo) 
respecto a la ecuación de la FAO Penman-Monteith (PM). 


Materiales y métodos 


La zona arrocera al sur del Valle del Cauca se extiende aproximadamente 
entre los 3% 5' 0” a los 3% 20' 0” de latitud norte y entre los 76% 10' 0” y 
76% 50' 0” de longitud oeste. La elevación de la zona es de unos 1 000 
msnm; predomina una temperatura media de 23 *C, con una oscilación 
diaria de 11 *C; húmeda relativa del 73 al 80 %; brillo solar de 138 a 200 
horas mensuales (Valencia, García, 8 Montero, 2017). Según Kóppen, 
corresponde al clima tropical lluvioso de bosque (Martínez 8 Mendivelso, 
2004: p. 53) (Figura 1). 
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LEYENDA 


Dpto. Valle Del Cauca 
Microlocalización [| 

Estaciones Climatológicas e 
Curvas de Nivel ----- 


Barranquilla 
O Giz, 


Cartagena 


3 PS E 
390 3"13'30"N 3'18'0N 3"22'30"N 


34'30"N 


Figura 1. Localización del área de estudio. 
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Se usaron los registros diarios de la temperatura máxima, 
temperatura mínima, humedad relativa y brillo solar de cinco estaciones 


climáticas por un periodo de estudio de 1989 a 2016 (Tabla 1). 


Tabla 1. Descripción de las estaciones climáticas. 


Altitud . ñ 
Latitud (N) Longitud (W) 
(msnm) 
Ciat Quilichao PEA 03% 03“ 76% 30“ 
Ingenio Bengala 1 000 03 15' 76% 24" 


977 
La Independencia 03 11” 8” 76% 34" 9” 
Ingenio Miranda 1 128 03% 16" 76% 20" 


Modelos para estimar la evapotranspiración de 
referencia diaria (Eto) 


Se calculó la ETo con los cinco modelos reconocidos a nivel diario en el 


programa Excel. A continuación se enlistan los modelos. 
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Ecuación Penman Monteith (PM) 


El método Penman-Monteith es considerado el modelo estándar y el más 
preciso para estimar ETo (Allen et a/., 1998). Fue desarrollado con la 
definición del cultivo de referencia, que es cultivo hipotético de altura de 
0.12 m, con una resistencia superficial de 70 s m*1 y un albedo de 0.23, 
y que representa la evapotranspiración de una superficie extensa de pasto 
verde de altura uniforme, creciendo activa y adecuadamente regado 
(Allen, Pereira, Raes, € Smith, 2006): 


900 
0.408A (Rn—G) +1 77 (es ea) 
A+y(1+0,34u>) 


ET0pm = (1) 
EToprm: evapotranspiración de referencia (mm día"!). 

Rn: radiación neta en la superficie del cultivo (MJ m-? d 71). 

G: flujo de calor de suelo (MJ m-? d 71). 

T: temperatura media del aire a 2 m de altura (*C). 

uz: velocidad del viento a 2 m de altura (m s 71). 

Es-ea: déficit de presión de vapor (kPa). 


A: pendiente de la curva de presión de vapor a saturación (KPa* C 71), 
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y: constante psicrométrica (Kpa* C*1). 


El procesamiento de los datos por ser llevado a escala diaria el valor 
G = 0 para el cultivo de referencia es según Allen et a/. (2006). 


Ecuación Hargreaves-Samani (HS) 


El modelo de Hargreaves-Samani es un modelo más simple que requiere 
de variables meteorológicas, temperatura (media, máxima y mínima), 
radiación incidente (Xu 8 Singh, 2001) y radiación solar extraterrestre 
(Ra) para estimar la ETo para una latitud y día determinados; Ra puede 
obtenerse a partir de tablas o calcularse mediante un conjunto de 


ecuaciones que utilizan la temperatura : 


ETO ys = 0.0023(Tmedia + 17.8) Tmax = Tmin Ra (2) 


Ecuación de Priestley-Taylor (PT) 
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Es un modelo semiempírico que requiere de menos parámetros para el 
cálculo de la evapotranspiración, comparado con el modelo de Penman- 
Monteith (Cervantes, Arteaga, Vázquez, Ojeda, 8 Quevedo, 2013). Esta 
metodología involucra la radiación solar y manifiesta menos errores 
respecto del método estándar (Pereira, 2004). Los elementos necesarios 


para su cálculo son: 


A: pendiente de la curva de vapor de saturación (KPa* C 71), 
A: coeficiente empírico equivalente a 1.26. 

y: constante psicrométrica (KPa* C-1), 

R,,: radiación neta del cultivo (MJ m2 d -1). 


G: flujo de calor de suelo (MJ m-2 d 1), por ser a escala diaria el valor G 


= 0 para el cultivo de referencia según Allen et a/. (2006): 


ETO pr = a+ ||» (n-0) (3) 


Ecuación de Linacre (LN) 
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Este modelo se deriva de Penman, estimando la evapotranspiración a 
través de datos geográficos (latitud y altitud) y de temperatura (Carvalho, 
Rios, Miranda, 8 Castro, 2011): 


(500 Tm 


do Emo )+15:(T-TPR) 


0-7 (4) 
EToLN, en mm día"!, 

Tm =T + 0.006 * h. 

h: altura del lugar, en msnm. 

TPR: temperatura del punto de rocío, “C. 

A: latitud del lugar (grados decimales). 


T: temperatura media, “C. 


Ecuación de García y López (GL) 


El modelo propuesto se adapta a las latitudes comprendidas desde los 150 
N hasta los 150 S. Las variables utilizadas son la temperatura media en 
OC y la humedad relativa en % (García 8 López, 1970): 
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7.45*T 
EToGL = [121 x 10 234.7+7 x (1— 0.01 x HR)|+0.21x T—2.3 (5) 


Modelo de regresión lineal 


Se utilizaron cinco estaciones para definir los diagramas de dispersión y 
la relación que existe entre los datos calculados con el modelo de PM y 
los estimados con las ecuaciones de HS, PT, LN y GL, de las cuales se 
obtienen los siguientes modelos lineales: 


ETOpm-us = 4 +bx*EToHS (6) 
ETOpy-pr = 4 + bx EToPT (7) 
ETOpm-.yN = 4+bx*EToOLN (8) 
ETOpy-c, = 4 +b*ETOGL (9) 


Índices estadísticos 
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Los índices utilizados para evaluar la bondad de ajuste de los modelos 


son: 


a) 


b) 


c) 


El coeficiente de determinación (R2) superior a 0.60 indica el grado en 
que la regresión explica la suma del total al cuadrado (Tagliaferre et 
al., 2012). 


La raíz cuadrada del cuadrado medio del error compara un valor 
predicho y un valor observado o conocido (RMSE) (Caí, Liu, Lei, Y 
Santos, 2007). 


Eficiencia de Nash-Sutcliffe (NSE): este coeficiente establece la 
relación entre el comportamiento real y el modelado. El dominio 
matemático del valor del coeficiente NSE es de —o a 1, donde la unidad 
representa la simulación perfecta. Un ajuste perfecto quiere decir que 
la varianza de los errores es cero; vale cero cuando la varianza de los 
errores es iguala la varianza observada, lo cual significa que el modelo 
produce estimaciones del promedio de las observaciones en todos los 
intervalos. Los valores negativos indican un desempeño malo del 
modelo (Magana-Hernández, Ba, Guerra-Cobian, €: Víctor, 2013). 


Diagrama de cajas 
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Se emplearon estos gráficos para identificar la forma en que se 
distribuyen la media y desviación estándar, y la manera en que se 


distribuyen los datos calculados portodos los modelos para cada estación. 


Resultados 


Se presenta el comportamiento de la evapotranspiración de referencia de 
algunos años hidrológicos representativos para la zona productora 
evaluada. El año 1999 se eligió como el más lluvioso (10 513 mm año“! 
para la zona), 2012 fue un año intermedio (7 791 mm año?) y 1992 un 
año seco (5 544 mm año”). 


De hecho, PT y GL presentan resultados similares en el año lluvioso; 
pero en el año intermedio y seco, GL presentó valores superiores a PT, 
pero inferiores a HS. Otra observación considerable es que el método de 
LN mostró estimaciones superiores en los meses de julio, agosto y 
septiembre, ya que la temperatura es superior en esta época del año e 


influye en la estimación por este método (Figura 2, Figura 3, Figura 4). 
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Figura 2. ETo promedio mensual para un año lluvioso en la zona 


productora de arroz. 
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Figura 3. ETo promedio mensual para un año intermedio en la zona 


productora de arroz. 
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Figura 4. ETo promedio mensual para un año seco en la zona 


productora de arroz. 


Se observó que para las estaciones Ing. Bengala y Ciat Quilichao se 
obtuvieron estimaciones superiores a 0.5 en el coeficiente de 
determinación en todos los modelos. De hecho, las estimaciones entre los 
modelos Penman-Monteith y Priestley-Taylor presentó valores de 0.597 
de R2 para Ing. Bengala y 0.834 para Ciat Quilichao (Figura 5) y valores 
inferiores a 0.35 de R2 para Ing. Cauca, Independencia e Ing. Miranda. 
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Ciat Quilichao 


y = 0.9413x+ 0.3996 


ETo Priestley-Taylor (mm) 


0 1 2 3 4 9 
ETo Penman-Monteith (mm) 


Figura 5. Regresión Lineal Entre PM y PT para la estación Ciat 
Quilichao. 


El modelo HS exhibió coeficientes de determinación bajos de la 
siguiente forma: valores de R2 ¡igual a 0.376, 0.290 y 0.268 para las 
estaciones Ing. Cauca, Independencia e Ing. Miranda, respectivamente, 
no así para la estación Ing. Bengala, donde su coeficiente de 


determinación fue de 0.8264 (Figura 6). Asimismo, el modelo de García y 
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López presentó un rango de valores de R? entre 0.42 y 0.29 para las tres 


estaciones mencionadas antes. 


Ingenio Bengala 


9 
E 8 
E 
al 
EG 
E 
hs 
174) 
94 
49] 
9 3 
Bo 
5 2 o y = 0.9047x+0.9278 
o; R?= 0.8264 
E 
uu 

0 

1.2 2.2 3.2 4.2 5.2 6.2 7.2 8.2 9.2 


Eto Penman-Monteith (mm) 


Figura 6. Regresión lineal entre PM y HS para la estación Ing. Bengala. 


Linacre se manifestó con los valores más bajos de R2 para 
Independencia e Ing. Cauca (Figura 7), aunque para Ing. Miranda exhibió 
un valor de 0.330 de R2. 
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ETo Penman-Monteith (mm) 
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Los intervalos de confianza en el que se estimaron los coeficientes 


de las regresiones lineales están basados en un 95 % de confianza (Figura 


5, Figura 6 y Figura 7, Tabla 2), y se observa que los valores del intercepto 


“a” y “b” están dentro de los límites permitidos para dicho nivel de 


confianza. 
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Tabla 2. Intervalo de confianza al 95 %. 
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Intervalo de Intervalo de 
confianza a confianza b 
intercepto| Pendiente 
Modelos| Estación Límite | Límite | Límite | Límite 
inferior| superior | Inferior | superior 
a a b b 


0.3996 0.9413x| 0.371 0.429 0.933 0.949 


0.9278 | 0.9047x | 0.892 0.964 0.897 0.913 
E ala 
LN-PM Az 2.8640 0.1578x | 2.829 2.899 0.148 0.168 


HS: ecuación de Hargreaves-Samani 


PT: ecuación de Priestley-Taylor 
LN: ecuación Linacre 


PM: Penman Monteith 


Para este estudio, el error promedio del modelo HS fue de 1.38 mm 
día"!, lo cual se exhibe como aceptable; aunque se evidencia como el error 
promedio más alto de todos los modelos, respecto a PM (Tabla 3). A su 
vez, este índice permite presentar que el modelo de Linacre puede ser un 
sustituto del método PT, pues su RMSE promedio de las cinco estaciones 
para el área de estudio fue 0.71 mm día"! (Tabla 3). 
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Tabla 3. Raíz cuadrada del cuadrado medio del error (RMSE) de la ETo- 
PM respecto a HS, PT, GL y LN para cada estación 


PM-Ing. | PM-Ing.| PM-Ciat 
Bengala | Cauca |Quilichao | Independencia | Miranda 


PT 0,35 0.59 0.36 1.13 0.76 


AA 


Ml ill EA 


En cuanto a la eficiencia de Nash-Sutcliffe (NSE), se halló poca 
precisión en la estimación de la ETo en todos los modelos por estaciones. 
No obstante, se identifican dos resultados positivos significativos: PM-PT 
para Ciat Quilichao y PM-LN para Ing. Bengala (Tabla 4); no hay duda de 
que el modelo PT en la estación Ciat Quilichao es el único modelo que 
obtuvo un valor de eficiencia positivamente alto y tuvo menor error (Tabla 
3 y Tabla 4). 


269 
Tecnología y ciencias del agua, ISSN 2007-2422, 13(1), 246-291. DOI: 10.24850/j-tyca-2022-01-06 


Ed 2022, Instituto Mexicano de Tecnología del Agua 
Tecnología y 


Ciencias¿Agua Open Access bajo la licencia CC BY-NC-SA 4.0 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


Tabla 4. Eficiencia de Nash-Sutcliffe (NSE) de la ETo-PM respecto a HS, 
PT, GL y LN para cada estación 


PM-Ing. | PM-Ing. | PM-Ciat 
Bengala | Cauca | Quilichao | Independencia | Miranda 


ld o MN dl 


PT -0.70 -2.06 0.74 -3,84 “2.97 


Diagrama de cajas 


Los resultados obtenidos de los análisis anteriores lograron identificar que 
el modelo Priestley-Taylores uno de los métodos que tiene alta afinidad 
con el método estándar; pero otra característica relevante de la discusión 
es la distribución del conjunto de datos de cada modelo por estación 
(Figura 8). 
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8 
o 
7 
sé Ml HS: Hargreaves-Samani 
E ; Ml PT: Priestley-Taylor 
o Mi GL: García-López 
mu 4 
, El LN: Linacre 


El PM: Penman Monteith 


Modelos ETo 


Figura 8. Estimaciones de la ETo para la estación Ciat Quilichao. 


Existe una gran variabilidad en las estimaciones de la ETo diaria, 
por lo que es notorio la posición desigual de las cajas entre sí, y tras las 
estimaciones de la NSE se identificó que para la estación independencia 
no se describe cuantitativamente la precisión de los modelos usados 


contra el modelo estándar (Figura 9). 
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Modelos ETo 


Figura 9. Estimaciones de la ETo para la estación Independencia. 


De la estación Ing. Miranda, los modelos LN y PT tienen una 
varianza de 0.12 y un comportamiento similar en su desviación estándar 
de 0.34 y 0.35 mm día”?*, con un conjunto de datos que coinciden con las 
estimaciones hechas por el método estándar (Figura 10). De hecho, el 
modelo LN-PM presentó mejores resultados que los demás modelos para 


esta estación (Tabla 3). 


272 
Tecnología y ciencias del agua, ISSN 2007-2422, 13(1), 246-291. DOI: 10.24850/j-tyca-2022-01-06 


2022, Instituto Mexicano de Tecnología del Agua 


Tecnología y 
Ciencias¿Agua Open Access bajo la licencia CC BY-NC-SA 4.0 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 
ING.MIRANDA 
6.5 
6 j o - PS 
5.5 5 8 $ 
8 
; j 


Ml HS: Hargreaves-Samani 
ll PT: Priestley-Taylor 


GL: García-López 
El LN: Linacre 
El PM: Penman Monteith 


ETo (mm) 
Dn a 


uu 
o 


Modelo ETo 


Figura 10. Estimaciones de la ETo para la estación Ing. Miranda. 


La estación Ing. Bengala exhibe que las estimaciones hechas por el 
modelo PM y HS tienen una varianza de 0.281 y 0.278, respectivamente, 
y una desviación estándar similar de 0.53 mm día”? (Figura 11). Sin 
embargo, para este caso fue el mismo método estándar el que presentó 
más observaciones que se alejan de manera poco usual del resto de los 
datos. En cambio, los modelos PT, GL y LN en su conjunto de datos están 
más cercanos a la media; a su vez, obtuvieron los errores más bajos para 
esta estación con el índice RMSE. 
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Figura 11. Estimaciones de la ETo para la estación Bengala. 


Los valores obtenidos de la raíz cuadrada del cuadrado medio del 
error para la estación Ing. Cauca fueron inferiores a 1 mm para los 
modelos PT, GL y LN (Figura 12, Tabla 3), los cuales concuerdan con el 
95 % de las estimaciones realizadas con el modelo GL con el modelo 
estándar para la estación Ing. Cauca. En cambio, el 75 % de los resultados 
obtenidos del modelo de Linacre coindicen con el 50 % de los valores 
obtenidos del modelo PM. Del mismo modo, Hargreaves-Samani posee la 
media más alta con valores atípicos y extremos, que comparando éste 
con el modelo estándar se podría indicar que HS no es un método 
conveniente para dicha estación, ya que PM los considera dentro de su 


varianza. 
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Figura 12. Estimaciones de la ETo para la estación Ing. Cauca. 


Los gráficos exhibidos anteriormente permiten identificar que los 
modelos más cercanos al método estándar para las cinco estaciones 
fueron el PT y LN, pues presentaron estadísticas idénticas, mejores 


ajustes o menores errores. 


Discusión 
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El modelo PM a nivel general se manifestó en un intervalo de 3 a 3.8 mm 
mes”1*y HS se exhibió como el método que presenta promedios mensuales 
superiores en todos los años elegidos. Enel año lluvioso, junio y diciembre 
se presenta como los meses con los promedios más bajos en todos los 
métodos; esto acontece por el fenómeno de la zona de convergencia 
intertropical, que influye en el comportamiento bimodal de la precipitación 
para el país (Figura 2). 


Marzo y agosto del año seco presentan los promedios más altos y 
entre métodos los valores tienden a ser similares (Figura 4); de acuerdo 
con el registro del IDEAM, agosto es el mes que mayor evaporación 
registra en los datos históricos. 


Un estudio hecho por Gunston y Batchelor (1983) para comparar 
los modelos PM y PT estimó la evapotranspiración de referencia a nivel 
mensual en los países tropicales, y registró un desempeño aceptable de 
R2 = 0.87 para 126 meses húmedos; R2=0.90 para 86 meses de lluvias 
intermedias, y R2?2=0.48 en 148 meses de sequía. Otro ejemplo fue 
ejecutado por Qiu et al. (2019), en la ciudad Nanjing, de la provincia 
Jiangsu en China, cuyo modelo PT estimó razonablemente la ETo para la 
zona de producción de arroz con unos coeficientes de determinación de 
0.92 y 0.96 para dos estaciones. Este modelo ha sido exitoso, después 
del método PM, gracias a la radiación neta (Rn), que es una variable 


meteorológica que afecta a la ET-referencia a corto plazo en el sistema 
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de producción de arroz. No obstante, cuando una estación meteorológica 
tiene un conjunto de datos limitados de temperatura y radiación solar, la 
mejor opción para ese caso parece ser el uso de PT, que es una 
simplificación del método original Penman-Monteith (Senthelas, Gillespie, 
8: Santos, 2010). 


Respecto al modelo HS, Awal, Habibi, Fares y Deb (2020) indican 
que la Elo diaria estimada con el modelo Hargreaves-Samani en 
diferentes estaciones del oeste de Texas, EUA, reveló más del 0.7 de 
coeficiente de determinación entre este modelo respecto al de PM; es 
decir, que existe un alto nivel de concordancia entre ambos métodos para 


esta zona. 


En cambio, los resultados del modelo GL, a pesar de ser bajos son 
considerados aceptables, pues este modelo ha sido empleado en otros 
estudios a nivel mensual y anual (Sánchez, 2000); es una opción para 


cuando no se pueden emplear otros métodos debido a la carencia de 


registros de brillo solar. 


En este sentido, los resultados del modelo LN se corroboran con los 
encontrados por Issaka et al. (2017), quienes indicaron que este modelo 
se posicionó en último lugar con un rendimiento moderado de R2 = 0.47 
entre la evaluación de seis modelos respecto al estándar, cuyos datos 
diarios utilizados fueron los de la estación meteorológica del Aeropuerto 
Internacional de Doha, Qatar, en un período de 30 años. Otro ejemplo es 
el de Araújo, Oliveira, Araújo, Ledo y Silva (2010), quienes evaluaron el 
desempeño de los métodos de estimación de la evapotranspiración de 
referencia en las ciudades del estado de Ceará, encontrando un 
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desempeño aceptable para este método de R? = 0.934. El rendimiento de 
este modelo se debe a que utiliza sólo datos de temperatura del aire y 
altitud como la variable del sitio, limitando la representatividad de las 
condiciones climáticas para propósitos de estimación de la ETo (Pereira, 
Yanagi, Mello, Silva, 8. Silva, 2009). 


De los 20 modelos lineales desarrollados, se identificó que el valor 
de la pendiente es estadísticamente diferente de cero, pero no mayor a 
uno, y por lo tanto la regresión es significativa. De hecho, el valor “p” fue 
similar(< 2e-16) para todas las regresiones; es un valor inferior al valor 


de “p” crítico (0.05). 


De los valores obtenidos del índice RMSE entre los modelos PM-HS 
se identifican similitudes expresadas por Amatya, Skaggs y Gregory 
(1995), quienes obtuvieron un intervalo de errores de RMSE de 0.87 a 
1.37 mm día"1; en cambio, el error promedio de RMSE entre el modelo HS 
y PM para 12 estaciones en Khuzestan, suroeste de Irán, fue de 1.17 mm 
día"! (Rahimi, 2008). Sin embargo, Lavado, Lhomme, Labat, Loup-Guyot 
y Boulet (2015) presentan 0.50 mm día"? como error de la media del 
modelo HS respecto al modelo de referencia PM para ocho estaciones en 
Perú. Si bien es notorio, la estación Independencia exhibió los resultados 
más altos en este índice en los modelos HS, PT y GL contra PM; esto 
permite enunciar que tales métodos se pueden recomendar para el cálculo 
de la ETo en otro tipo de escala para dicha estación (Martínez, Boueri, € 
Escalona, 2005). 


Los resultados de LN son valores similares a los reportados por 


Lavado et al. (2015), quienes exhiben valores de 0.42 a 1.05 mm día"! 
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para ocho estaciones en la cuenca amazónica andina. Esto ha sido posible 
gracias a que se requiere sólo la temperatura media (T) y el punto de 
rocío (Td), y es un modelo útil en una amplia gama de climas expuesto 
por Linacre (1977); no es necesario medir el punto de rocío si se dispone 
de temperaturas extremas diarias, pues éstas permiten hacer una 


estimación aproximada. 


Para el alcance de este trabajo, el valor obtenido de la NSE para 
Ciat Quilichao resulta ser mucho mayor que el reportado por Celestin, Qi, 
Li, Yu y Cheng (2020), quienes evaluaron 32 ecuaciones de 
evapotranspiración de referencia contra el método Penman-Monteith en 
el corredor Hexi, en el noroeste de China; también concluyen que el 
modelo HS resulta ser un buen sustituto de PM a escala mensual y para 
dicho lugar. Por otro lado, del modelo Linacre no se halló trabajos que 


indaguen sobre la bondad del modelo. 


Diagrama de cajas 


El modelo PM exhibe un rango de 3 a 4 mm día"! de la ETo estimada para 
la estación Ciat Quilichao, la cual coinciden con el 50 % de los resultados 
obtenidos por el modelo PT. Ambos modelos exhiben correspondencia 
similar para esta estación, ya que los respectivos valores de varianza son 
de 0.50 mm2 para PM y 0.53 mm? para PT, y la desviación estándar es 
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de 0.71 mm día"! para PM y 0.73 mm día"! para PT, lo que permite inferir 
que en ausencia de datos para estimar este parámetro con el modelo de 
PM se puede emplear el PT, tal cual se corroboró con el índice de la 
eficiencia de Nash-Sutcliffe (Tabla 4 y Figura 8), similitud que también 
divulgan Lu, Sun, McNulty y Amatya (2005) al comparar seis modelos 
diarios para calcular la ETo respecto al modelo PM en el sureste de EUA. 


De hecho, Linacre se puede considerar como un tercer modelo para 
estimar la ETo en la estación Independencia, ya que ha manifestado un 
intervalo de resultados similares al del PM. En cambio, de los otros 
modelos se puede mencionar que no coinciden con PM, pues la varianza 
para GL, PT y HS es de 0.60, 0,40 y 0.53 mm, respectivamente (Figura 
9). De modo comparativo, Vicente-Serrano eta/. (2014) exhiben que para 
la época de otoño, el modelo Linacre mostró valores de ETo igual y más 
altos que el método PM, pero se observaron valores inferiores estimados 
de LN para la sección de verano respecto al método estándar. Esto es 
posible, ya que cada ecuación fue creada para unas condiciones climáticas 
específicas que pueden diferir de las presentadas por cada estación; así, 
el comportamiento de cada modelo es muy cambiante dependiendo del 
lugar y época del año en la que se utilice (Sánchez, 2000). 


No obstante, en la estación Ing. Miranda, el modelo PM presenta un 
conjunto de datos entre 1.7 y 5.1 mm día"*, con pocos valores extremos 
entre 5.2 y 5.5 mm día"!*; eso indica que el uso de este modelo para esta 
estación exhibe algunos valores alejados de su media en más de dos 
unidades (Figura 10). Una investigación similar realizada en Macedonia 
Occidental, Grecia, por Efthimiou, Alexandris, Karavitis y Mamassis 
(2013), presentó que los modelos LN y PT reportaron el 79 y 89 % de 
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coincidencia en las estimaciones diarias de 34 años de la ETo respecto a 
las del Penman-Monteith. 


Tanto el Modelo PM y Hs sirven para caracterizar dicha variable en 
la estación Ing. Bengala; esto se confirma por el coeficiente de 
determinación entre ambos métodos que fue de 0.826 (Figura 6). PM se 
ve influenciado por varios parámetros que, de acuerdo con los 
experimentos hechos por Tyagi, Sharma y Luthra (2000), existe una 
probabilidad del 93.4 % que un modelo diferente al de PM tenga 
coincidencia con él (Figura 11); pero Lang et al. (2017) aluden que la 
ecuación de HS se puede usar más fácilmente que PT para estimar la ETo 
en el suroeste de China, ya que la primera requiere menos parámetros 


para estimarla ETo que la segunda. 


En la estación Ing. Cauca, los modelos PT y GL se caracterizan por 
presentar 0.43 y 0.44 mm día"? de desviación estándar, pero con una 
media y conjunto de datos diferentes. Gong, Xu, Chen, Halldin y Chen 
(2006) exponen que las variables climáticas en las mediciones diarias de 
la ETo entre el modelo PM y GL manifiestan que la humedad relativa es el 
factor meteorológico más sensible, seguido de la velocidad del viento, la 
radiación solar y la temperatura. En la Figura 12, el método que es más 
parecido al de PM es el de LN. El modelo HS presenta la mayor dispersión 
en la distribución de los datos; casos similares los reportan Pérez y 
Castellví (2002), quienes enuncian que este modelo exhibe los mayores 


valores atípicos en el cálculo de la ETo en la región de Cataluña en España. 


En síntesis, la relación de las ecuaciones Penman-Monteith y 


Priestley-Taylor expone su validez cuando se disponen de registros 
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necesarios para su aplicación a nivel diario. Dos anotaciones importantes 
las expresan Vicente-Serrano et al. (2014), quienes indican que los 
modelos basados en las temperaturas como Hargreaves-Samani suelen 
subestimar la ETo para varias estaciones de un solo lugar y que algunos 
métodos pueden mostrar la mejor concordancia con el modelo PM, pero 
podría tener la menor concordancia en una estación diferente. Sin 
embargo, es importante manifestar que los modelos restantes 
constituyen una herramienta admisible cuando no se cuenta con los datos 
climatológicos suficientes o de procedencia confiable, tal y como lo 
señalan Allen et al. (1998), quienes señalan que se pueden sugerir 
diversos métodos para el cálculo de la ETo con base en el tipo de 


información climática que se tenga. 


Conclusiones 


El modelo Penman-Monteith de la FAO estima la evapotranspiración de 
referencia diaria (ETo) al considerar un conjunto completo de datos 
meteorológicos y, a su vez, es la principal restricción para su uso en 
lugares donde faltan datos; situación que se da en el sur del 


departamento del Valle del Cauca, y con pocas estaciones meteorológicas, 
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donde es frecuente la falta de datos sobre velocidad del viento, la 
humedad relativa y la radiación solar. No existe ningún estudio previo que 
evalúe el rendimiento de los diversos métodos para estimar la 


evapotranspiración de referencia diaria en esta región de Colombia. 


El modelo Priestley-Taylores un buen método entre los estudiados; 
manifestó una alta similitud con los valores obtenidos respecto al modelo 
estándar Penman-Monteith con el índice de la raíz cuadrada del cuadrado 
medio del error (RMSE) para todas las estaciones. Con el índice de la 
Eficiencia de Nash-Sutcliffe sólo presentó un valor altamente aceptable 
para la estación Ciat Quilichao. Otro modelo útil para estimar la 
evapotranspiración de referencia diaria es el de Linacre; al igual que 
Priestley-Taylor, ambos métodos son ideales, ya que son modificaciones 
del modelo estándar. 


Existe evidencia de que el uso a nivel diario del modelo García y 
López es limitado, a pesar de ser una ecuación adaptada para el trópico, 
cuyos factores climatológicos son de fácil adquisición. Respecto al modelo 
Hargreaves-Samani, se presenta como un método aceptable para la 
estación Ing. Bengala; sin embargo, este modelo junto al de García y 
López no son ideales para la estación independencia, ya que los resultados 
obtenidos por el método RMSE Y NSE no son acordes con lo esperado en 
los alcances de este trabajo. 
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